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Abstract. We present a large bank of chemical profiles and pulsation periods suited 
for asteroseismological studies of ZZ Ceti (or DAV) variable stars. Our background 
equilibrium DA white dwarf models are the result of fully evolutionary computations 
that take into account the complete history of the progenitor stars from the ZAMS. 
^3 ■ The models are characterized by self-consistent chemical structures from the centre 

to the surface, and cover a wide range of stellar masses, effective temperatures, and 
H envelope thicknesses. We present dipole and quadrupole pulsation g-mode periods 
comfortably covering the interval of periods observed in ZZ Ceti stars. 

Complete tabulations of chemical profiles and pulsation periods to be used in as- 
teroseismological period fits, as well as other quantities of interest, can be freely down- 
loaded from our website (http : //www. fcaglp .unlp . edu. ar/evolgroup). 



1. Introduction 

White dwarf asteroseismology is a powerful astrophysical tool that fully exploits the 
comparison between the observed pulsation periods in white dwarfs and the periods 



, computed for appropriate theoretical models. It allows us to infer details of the origin, 

internal structure and evolution of white dwarfs (Winget & Kepler 2008; Althaus et 
al. 2010a). ZZ Ceti (or DAV) stars constitute the most numerous group of degenerate 
variable stars. They are otherwise normal DA (H-rich atmospheres) white dwarfs that 
exhibit g(gravity)-mode pulsations. Recently, our group La Plata Stellar Evolution 
and Pulsation Research Group has performed for the first time a detailed asteroseis- 
mological study on an ensemble of 44 bright ZZ Ceti stars by employing fully evo- 
lutionary (that is, non static) DA white dwarf models (Romero et al. 2012, Romero 
2012). Our asteroseismological approach basically consists in the employment of a 
large suite of detailed stellar models characterized by very accurate and updated physi- 
cal ingredients. These models were produced by computing the complete evolution of 
the progenitor stars. We have applied successfully this approach to the hot GW Vir (or 
DOV) stars in the past (see Corsico et al. 2007a, 2007b, 2008, 2009). Since the final 
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chemical stratification of white dwarfs is fixed in prior stages of their evolution, the 
evolutionary history of progenitor stars is of utmost importance in the context of white 
dwarf asteroseismology. 

Here, we present a complete set of pulsational results that can be useful to perform 
asteroseismological studies of ZZ Ceti stars. These include the internal chemical pro- 
files and the run of the Brunt- Vaisala and Lamb frequencies, as well as a large set of 
adiabatic pulsation periods, time-averaged oscillation kinetic energies, and first order 
rotational splitting coefficients. All these quantities can be freely downloaded from our 
website. 

This database can be used in two ways: (1) carrying out period fits using directly 
the bank of periods we provide, corresponding to our set of equilibrium DA white dwarf 
models, or alternatively, (2) by scaling our internal chemical profiles to the structure of 
other independent DA white dwarf models, and then by performing period fits using 
the pulsation periods computed for such models. 

Below, we summarily describe the input physics of our models and the pulsation 
computations, and then we present the format of our database. 

2. Input physics and evolutionary computations 

The state-of-the-art DA white dwarfs evolutionary models employed have been com- 
puted with the LPCODE evolutionary code. Details of this code and the input physics 
that this code includes can be found in Althaus et al. (2010b), Renedo et al. (2010) and 
references therein. Below, we briefly enumerate the main physical ingredients: 

• For the high-density regime characteristics of white dwarfs, we have used the 
equation of state (EoS) of Segretain et al. (1994). 

• For the low-density regime, we used an updated version of the EoS of Magni & 
Mazzitelli (1979). 

• Radiative opacities are from the OPAL project (Iglesias & Rogers 1996), in- 
cluding carbon- and oxygen-rich composition, supplemented at low temperatures 
with the molecular opacities of Alexander & Ferguson (1994). 

• Conductive opacities are those of Cassisi et al. (2007). 

• Neutrino emission rates are taken from Itoh et al. (1996) and Haft et al. (1994). 

• The 12 C(a, y) 16 reaction rate, of special relevance for the C-0 stratification of 
the white dwarf, is that of Angulo et al. (1999). 

• Convection has been modeled with the formalism of the mixing-length theory as 
given by the ML2 parametrization (Tassoul et al. 1990). 

• Treatment of the chemical profiles. We have considered the distinct physical 
processes that are responsible for changes in the chemical abundance distribution 
during white dwarf evolution: 

- Element diffusion: it strongly modifies the chemical composition profile 
throughout their outer layers. We have considered gravitational settling as 
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well as thermal and chemical diffusion of ! H, 3 He, 4 He, 12 C, 13 C, 14 N, and 
16 0. The treatment of time-dependent diffusion is based on the multicom- 
ponent gas treatment presented in Burgers (1969). Diffusion becomes oper- 
ative once the wind limit is reached at high effective temperatures (Unglaub 
& Bues 2000). As shown in previous studies, diffusion in the white dwarf 
envelope is a key ingredient as far as the mode-trapping properties of ZZ 
Ceti stars is concerned. 

- Abundance changes resulting from residual nuclear burning (mostly dur- 
ing the hot stages of white dwarf evolution) have been taken into account. 
Nuclear burning fixes the maximum value of Mh that is expected in a DA 
white dwarf. 

- Chemical re-homogenization of the inner carbon-oxygen profile induced by 
Rayleigh-Taylor instabilities has been considered following Salaris et al. 
(1997). These instabilities arise because of the positive molecular weight 
gradients that remain above the flat chemical profile left by convection dur- 
ing core helium burning. 

3. Pulsation computations 

We have employed the pulsation code described in Corsico & Althaus (2006), which is 
coupled to the LPCODE evolutionary code. The pulsation code is based on a general 
Newton-Raphson technique that solves the full fourth-order set of equations governing 
linear, adiabatic, nonradial stellar pulsations following the dimensionless formulation 
of Dziembowski (1971) (see Unno et al. 1989). The pulsation code provides the dimen- 
sionless eigenfrequency uj^ (k being the radial order of the mode) and eigenfunctions 
yi, • • • ,y4. From these basic quantities, the code computes pulsation periods (I~U), time- 
averaged oscillation kinetic energies (£"&)> rotation splitting coefficients (Ct), weight 
functions (Wit), and variational periods (IIJ) for each computed eigenmode (see Ap- 
pendix A of Corsico & Althaus 2006 for the definition of these quantities). Usually, the 
relative difference between IT^ and II* is lower than ~ 10 -4 , which represents the typi- 
cal uncertainties in our theoretical periods. The prescription we follow to assess the run 
of the Brunt-Vaisala' frequency (AO for a degenerate environment typical of the deep 
interior of a white dwarf is the so-called "Ledoux Modified" treatment (Tassoul et al. 
1990), appropriately generalized to include the effects of having three nuclear species 
(oxygen, carbon, and helium) varying in abundance in the same place of the star (a 
triple chemical interface). Fig. [T] displays an example of the chemical profiles of our 
models with a fixed stellar mass and effective temperature and for different thicknesses 
of the H envelope (upper panel), and the run of the Brunt-Vaisala and Lamb (for £ - 1) 
frequencies (lower panel). Details can be found in Romero et al. (2012) and Romero 
(2012). 

We present results for DA white dwarf models with stellar masses, effective tem- 
peratures and H envelope thicknesses in the ranges 0.525 < M*/M G < 0.878, 9000 < 
r eff < 14000 K, and -3.62 < log(M H /M H .) < -9.34, respectively. We have performed 
pulsation calculations on about (11 x 7 x 200) = 15400 DA white dwarf models. In 
this account, we have considered the number of stellar mass values (11), the number 
of thicknesses of the H envelope for each sequence (~ 7), and the number of models 
(« 200) with effective temperature in the interval 14 000 - 9000 K, respectively. For 
each model, adiabatic pulsation g-modes with I = 1 and 2 and periods in the range 
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Figure 1. Upper panel: example of the internal chemical profiles for H, He, 
12 C, and 16 of our DA white dwarf models. The cases shown in the figure corre- 
spond to models with a stellar mass of M» = O.593M , an effective temperature of 
r eff ~ 12 350 K, and H envelope thicknesses of log(M H /M.) = -3.93 (black), -4.28 
(red), -4.85 (green), -5.34 (blue), -6.33 (magenta), -7.34 (brown), -8.33 (dark 
green), and -9.33 (violet). Lower panel: the logarithm of the squared Bmnt-Vaisala 
and Lamb frequencies corresponding to the cases depicted in the upper panel. The 
correspondence between each bump-like feature in N 2 and the chemical transition 
regions of the models is clearly visible. 



80 - 2000 s have been computed. This range of periods corresponds (on average) to 
1 < k < 50 for t = 1 and 1 < k < 90 for I = 2. So, more than ~ 2 x 10 6 adiabatic 
pulsation periods have been computed. An example of our pulsation results is depicted 
in Fig. |2j where we show the forward period spacing, the oscillation kinetic energy, and 
the first order rotational splitting coefficients for the same models analyzed in Fig. Q] 
Note the gradual changes experienced by these quantities as we vary the value of the H 
envelope thickness. 

In Table Q] we show the grid of evolutionary sequences included in our database. 
The values of the stellar mass of our set of DA white dwarf models are shown in the 
upper row. The masses of H corresponding to the different envelope thicknesses for 
each stellar mass are shown in rows 2 to 9. Row 2 corresponds to the maximum value 
of the thickness of the H envelope for each stellar mass according to our evolutionary 
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Forward period spacing logtkinetic energy) rotation splitting coefficient 




500 1000 1500 2000 500 1000 1500 2000 500 1000 1500 2000 



Period [s] Period [s] Period [s] 

Figure 2. Pulsation quantities for i = 1 modes corresponding to the same models 
analyzed in Fig. Q] In each panel, the curves have been arbitrarily shifted upwards for 
clarity, except the lowest one. Left panel: forward period spacing (AIL; = IL + i - IL) 
vs periods; middle panel: oscillation kinetic energy {log(Ek)) vs periods; right panel: 
first order rotation splitting coefficient (Ck) vs periods. 



6 



Romero, Corsico, Althaus, Miller Bertolami 



computations ("canonical envelopes"). Rows 3 to 9 correspond to H envelopes thinner 
than the canonical ones. 



Table 1. 

MJM Q 0.5249 0.5480 0.5701 0.5932 0.6096 0.6323 0.6598 0.7051 0.7670 0.8373 0.8779 

log(^) -3.62 -3.74 -3.82 -3.93 -4.02 -4.12 -4.25 -4.45 -4.70 -5.00 -5.07~ 

-4.27 -4.27 -4.28 -4.28 -4.45 -4.46 -4.59 

-4.85 -4.85 -4.84 -4.85 -4.85 -4.86 -4.87 -4.88 -4.91 

-5.35 -5.35 -5.34 -5.34 -5.35 -5.35 -5.35 -5.36 -5.37 -5.41 -5.40 

-6.33 -6.35 -6.33 -6.33 -6.34 -6.34 -6.35 -6.35 -6.35 -6.36 -6.39 

-7.34 -7.33 -7.34 -7.34 -7.33 -7.35 -7.33 -7.35 -7.34 -7.36 -7.38 

-8.33 -8.33 -8.31 -8.33 -8.33 -8.33 -8.33 -8.34 -8.33 -8.34 -8.37 

-9.25 -9.22 -9.33 -9.33 -9.25 -9.34 -9.33 -9.34 -9.33 -9.34 -9.29 



4. Format of the files to be downloaded 

The files of our database can be downloaded from our website: 



http : //www . f caglp . unlp . edu . ar/evolgroup/TRACKS/PULSATIONS/PULSATIONS_DA/pulsations_cocor 

Once in the site, you can see that the files are organized in three separate tables, which 
have exactly the same form as Table Q] above. For a given table, each element is as- 
sociated to a specific sequence characterized by the corresponding values of M* and 
log(MH/M*). In the first table we provide results of periods (II^) and the other quanti- 
ties (AI~U, log(2?fc) and Ct) for the harmonic degree I = 1, in terms of the radial order k. 
The same information, but for the case ( = 2, is provided in the second table. Finally, 
the third table includes the squared Brunt- Vaisala (N 2 ) and Lamb {L 2 ( ) frequencies, the 
Ledoux term (B), and the chemical abundances (by mass) of hydrogen (Xh), helium 
(Xne), carbon (Xq), and oxygen (Xo) in terms of the radial coordinate (r/R*) and the 
outer mass fraction coordinate, [- log(l - M r /M*)]. 

4.1. The first table: periods and other quantities for £ - 1 

Each element of the table is a hyperlink that leads to a given tar gzipped file. You can 
download a specific tar gzipped file by simply clicking on a given element of the table. 
For instance, if you click at element -3.62 in the first column of the table, you will 
download the file: |ava_® . 52 5_3 . 623 . tgz| This file corresponds to the sequence with 
M* - 0.525M G and log(M H /M*) = -3.623. After you detar the file, you should obtain 
a sequence of files: 

|ava(S>(S>l_(S> . 525_3 . 623 . outl 
|ava6>6>2_6> . 52 5_3 . 623 . outl 
|avaQQ3_Q.525_3.623Tout| 
|ava(S>(S>4_(S>.525_3.623Tout| 



Each of these files corresponds to a different, decreasing effective temperature. For in- 
stance, in Fig. [3] we show the appearance of the file avaOO 1_0 . 5 2 5_3 .623. out The 
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Figure 3. The aspect of the file ava001_0 . 525_3 . 623 . out corresponding to a 
DA white dwarf model with M« = 0.525M G , LJL e = 1.9822 x 1(T 2 , r eff = 17 703 
K, t = 11212.811 x 10 6 yr, and log(M H /M») = -3.623. 
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heading contains the value of the stellar mass in solar units flM_*/M_ sun), the stellar lu- 
minosity in solar units ( ]L_*/L_sun ), the effective temperature in Kelvin (Teff [K] ), 
the age in 10 6 yr units (Time [le6yr] ), and the asymptotic period spacing in sec- 



onds (DPa [s]| ), computed as in Tassoul et al. (1990). Below the heading, the fol- 
lowing quantities are listed: the harmonic degree (1), the radial order (k), the periods 
(Period), the forward period spacings (DP), the logarithm of the oscillation kinetic 
energies (log(Ek) ), and the first order rotation splitting coefficients flC_kll ). 

A very important quantity such as the rate of period change of a given pulsation 
mode with radial order k (FL. = dUk/dt), can be easily computed from the age (r) and 
the period (Yl/J of a given model (J) and the corresponding values of the previous model 
U - 1): 

du k _ nj - nf 1 

dt tJ - W" 1 

Finally, we also give the option of downloading the complete set of periods (for all 
stellar masses, H envelopes and effective temperatures) at once in the hyperlink below 
the table. 



4.2. The second table: periods and other quantities for I - 2 

All the above explanation holds also for the case of the second table, but in this case, 
the results correspond to I = 2 pulsation modes. 



4.3. The third table: critical frequencies and chemical profiles 

Now, if you click (for instance) at element -3.62 in the first column of the table, you 
will download the file: par_® . 52 5 _3 . 623 . tgz. Again, this file corresponds to the 
sequence with M* = 0.525M© and log(M H /M*) = -3.623. After you detar the file, you 
should obtain a sequence of files: 
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out| 


par©92_9 


525. 


.3 


623 


out| 


par©93_© 


525. 


.3 


623 


out| 


par®®4_® 


525. 


.3 


623 


outl 



Again, each of these files correspond to a different, decreasing effective temperature. 
But, at variance with the cases of the first and second tables, in this case the interval 
in T e tf between consecutive files is quite large (~ 500 K). In Fig. |4] we show the ap- 
pearance of the file par®®l_®. 525 _3 . 623 . out[ The parameters at the heading are, 
in this case, the stellar mass, the stellar luminosity, and the effective temperature. Be- 
low the heading, the following quantities are listed: the normalized radial coordinate, 
(r/R_sun), the outer mass fraction coordinate (-logq), the squared Brunt- Vaisala fre- 
quency (IN" 21) . the squared Lamb frequency (L"2), the Ledoux term (BLedoux), and 
finally the abundance by mass of hydrogen (HI), helium (He4), carbon (C12), and oxy- 
gen MB- 
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Figure 4. The aspect of the file par0ISl_(S . 5 2 5_3 . 623 . out corresponding to a 
DA white dwarf model with M„ = O.525M , L„/L = 7.22 x 1(T 3 , T eS = 13 933 K, 
andlog(M H /M„) = -3.623. 
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5. Final remarks 

• Additional results for other values of the H envelope thickness that are not in- 
cluded in this database, can be obtained by request to the authors. 

• If you use this database and publishes your results, please cite the paper: "To- 
ward ensemble asteroseismology of ZZ Ceti stars with fully evolutionary mod- 
els", Romero, A. D., Corsico, A. H., Althaus, L. G., Kepler, S. O., Castanheira, 
B. G., Miller Bertolami, M. M. 2012, MNRAS, 420, 1462 

• Any comment/criticism that helps to improve this database will be greatly appre- 
ciated! 
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